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The separation of a heavy quark and antiquark pair leads to the formation of a tube of flux, 
or "string", which should break in the presence of light quark-antiquark pairs. This expected 
zero-temperature phenomenon has proven elusive in simulations of lattice QCD. We study mixing 
between the string state and the two-meson decay channel in QCD with two flavors of dynamical 
sea quarks. We confirm that mixing is weak and find that it decreases at level crossing. While our 
study does not show direct effects of internal quark loops, our results, combined with unitarity, give 
clear confirmation of string breaking. 
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I. INTRODUCTION 

In the absence of dynamical sea quarks, the heavy quark-antiquark potential is known quite accurately from 
numerical simulations of lattice quantum chromodynamics The potential is traditionally determined from the 
Wilson-loop observable, which is proportional to exp[—V{R)T] at large T. At large separation R, the potential V{R) 
rises linearly, as expected in a confining theory. In the presence of dynamical sea quarks the potential is expected 
to level off at large R, signaling string breaking. Thus far, no SU(3) simulation at zero temperature with light sea 
quarks has found clear evidence in the Wilson-loop observable for string breaking , even out to R ~ 2 fm. 

The reason string breaking has not been seen using the traditional Wilson- loop observable is now clear [^-0 . The 
Wilson loop can be regarded as a hadron correlator with a source and sink state (F) consisting of a fixed heavy quark- 
antiquark pair and an associated flux tube. The correct lowest energy contribution to the Wilson-loop correlator at 
large R should be a state "M" consisting of two isolated heavy-light mesons. However, such a state with an extra 
light dynamical quark pair has poor overlap with the flux-tube state, so it is presumably revealed only after evolution 
to very large T. To hasten the emergence of the true ground state, it is necessary to enlarge the space of sources to 
include both F and at least one M state. 
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Drummond demonstrated string breaking in a strong-coupling, hopping parameter expansion with Wilson quarks 
10,^. A number of numerical studies of theories less computationally demanding than QCD, including nonabelian 
theories with scalar and adjoint matter fields, found string breaking |^|l^. One study claims to have found string 
breaking in the absence of dynamical sea quarks by doing a transfer matrix calculation . 

In a full SU(3) simulation, until recently, string breaking has only been observed at nonzero temperature (close to, 
but below the deconfinement crossover) based on the Polyakov loop observable, which evidently has much better 
overlap with the M state. Some of us reported a preliminary low-statistics result for staggered quarks in 1999 p5| , 
and, last year, Pennanen and Michael announced evidence for string breaking at zero temperature using Wilson-clover 
quarks and a novel technique for variance reduction in computing the light quark propagator [ p^ . Duncan, Eichten, 
and Thacker found hints of a flattening static potential using a truncated determinant algorithm jlj]. 

In this paper we demonstrate string breaking in an SU(3) simulation with two flavors of dynamical sea quarks. 
Our simulations are done in the staggered fermion scheme on an archive of 198 configurations of dimensions 20'^ x 24, 
generated with the conventional one-plaquette Wilson gauge action at Q/g^ = 5.415 and two flavors of conventional 
dynamical quarks of mass am = 0.0125. At this gauge coupling and bare quark mass the lattice spacing is approxi- 
mately 0.163 fm (based on a measurement of the Sommer parameter jisll ro extracted from Wilson loops) with a pi 
to rho mass ratio of 0.358. These parameters were selected to give a relatively light quark, making pair production 
energetically favorable, and a large lattice volume (about 3.3 fm on a side and 3.9 fm in temporal extent) to allow 
ample room for string breaking. 

Our computational methodology is described briefly in Sec. |l|. In Sec. Ill wc justify our fitting ansatz. Finally, 
in Sec. |^ we present our results and conclusions. Two Appendices describe our formalism for random sources, and 
review the transfer matrix formalism we employ in our analysis. 



II. COMPUTATIONAL METHODOLOGY 



Our conventional Wilson loop is computed with APE smearing |19| of the space-like gauge links. Specifically, we 
used 10 iterations, combining the direct link with a factor 1 — a (in our case, a = 0.294) and six staples with factor 
a/6 with SU(3) projection after each iteration. In hamiltonian language the expectation value of this operator is the 
correlator Gff{RtT) between an initial and final state consisting of a static quark-antiquark pair separated by a 
fat string of color flux. Most of our results are obtained from on-axis Wilson loops with R ranging from 1 to 10, but 
we have two off-axis points at displacement (2,2,0) and (4,4,0) (plus permutations and reflections). Including other 
off-axis displacements might have been statistically useful 

We enlarge the source and sink space by including a meson-antimeson state M with an extra light quark located near 
the static antiquark and an extra light antiquark, near the static quark. To be precise, this state is the tensor product 
of a static-light meson operator and a static-light antimeson operator. As discussed in Appendix B, staggered flavor 
considerations make other choices more desirable close to the continuum limit, e.g. the tensor product of creation 
operators for a flux-tube state and a sigma meson. Our static-light meson construction makes the numerical analysis 
tractable and is adequate for studying mixing on coarse lattices. 

We use an extended source for the light quark in the static-light meson. The heavy quark position, on the other 
hand, is fixed and used to define the separation R. Specifically, the gauge-invariant source wavefunction at a site has 
support only on the site itself and on the second on-axis neighbors in all six spatial directions, connected to the central 
site by a product of the APE smeared links along the paths. Rather arbitrarily, the central site is given weight 2 
while the satellite sites have weight 1. Thus we also compute the additional correlation matrix elements Gmm{R,T), 
Gmf{R,T) and Gfm{R,T). They are diagramed in Figs. 0and||. 



FIG. 1. The static-light meson-antimeson pair contribution to the full QCD propagator. The wiggly lines denote the light 
quark propagator. Shown are the 'direct' and 'exchange' terms respectively. 
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FIG. 2. The string-meson correlation matrix element Gfm (and its hermitian conjugate Gmf)- The wiggiy line again denotes 
the light quark propagator. 

For the light quarks in the static-light mesons we use the same parameters as for the dynamical sea quarks. To 
reduce variance we generated "all-to-all" propagators for the light quark, using a Gaussian random source method 
(See Appendix Results reported here are based on 128 such sources per gauge configuration, which gave 
satisfactory statistics for our lattice volume. With this number of sources the variance due to fluctuations in random 
source was comparable to that due to fluctuations in the ensemble of gauge configurations. 

We analyze our correlators using an extension of the transfer matrix formalism of Sharatchandra, Thun, and 
Weisz [0 for staggered fermions, described in Appendix With our choice of local meson operators, discrete lattice 
symmetries, also discussed in Appendix B, require that all products of gauge links in the observables be assigned phases 
consistent with being viewed as the paths of heavy staggered fermions. For example, for the Wilson-loop operator, 
a hopping parameter expansion around an on-axis R x T rectangular path gives, in addition to the conventional 
Wilson-loop gauge-hnk product, a net phase factor (— independent of the staggered fermion Dirac phase 
conventions. (Included is a factor —1 for a single closed fermion loop.) This phase then controls the sign of the transfer 
matrix eigenvalue associated with the flux tube state. We use a similar construction to get the phases for the nonclosed 
gauge-link products in the diagrams of Figures |l| and ||. In all cases the Dirac phase convention for the gauge-link 
products must be consistent with that of the light quark. A consequence of this construction is that the correlation 
matrix elements involving off-axis gauge- link products must vanish when summed over symmetry-equivalent paths for 
net displacements R that have more than one odd Cartesian component. 

III. STAGGERED FERMION PAIR PRODUCTION 
A. Transfer matrix eigenvalues 

The heavy quark potential is defined as the ground state energy of the QCD hamiltonian with a static quark and 
antiquark separated by distance R. Operationally, we extract the ground state energy by fitting the time dependence 
of the correlation matrix elements in the same manner as one obtains hadron masses. Fundamentally, the potential 
is determined by the eigenvalues of the transfer matrix. To justify our fitting ansatz, therefore, we start from an 
analysis of the transfer matrix in the staggered fermion scheme. 

Transforming to temporal axial gauge and making a suitable choice of fermion phases, Sharatchandra, Thun, and 
Weisz showed that the staggered fermion transfer matrix is hermitian but not positive |po| . Then, it is convenient 
to use the eigenvectors of the transfer matrix as a basis for representing the correlation matrix. In terms of the 
(possibly negative) eigenvalues A„(i?) of the transfer matrix, our correlation matrix can therefore be written in 
spectral decomposition as 

N 

Gab{R,T) = ^Zi,(i?)ZB.(i?)[A,(i?)]^+i , (1) 

i=l 

where A, B refer to the flux tube F or meson-meson M states. The T -I- 1 power is natural, as we show in Appendix 
^. This result forms the basis for our fitting ansatz. To apply this decomposition to our results, it is essential, as we 
have done, that we treat the heavy quark lines as static staggered quarks, with all fermion phases included, and that 
source and sink operators are equivalent. 
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B. Static light propagator 



We first examine the single static-light meson correlator, shown in Fig. We find good fits to two spectral 
components, one with no phase oscillation in T, corresponding to a S'-wave light quark and a positive transfer- 
matrix eigenvalue, and the other, higher in energy, and with oscillating phase in T, corresponding to a P-wave light 
quark and a negative transfer-matrix eigenvalue. Fitting to a single nonoscillating exponential plus a single oscillating 
exponential over the range t = [2, 9] gives energies (defined as usual as log(|A|)) aEs = 0.7884(12) and aEp = 1.022(6) 
with /df = 2.7/4. The P-wave amplitude is suppressed by a factor of about 0.5 relative to the S'-wave amplitude. 




5 10 15 

T 

FIG. 3. Static light propagator with a nonoscillating S'-wave and oscillating P-wave component. The solid line connects the 
best fit values. 



C. Fitting form for the correlation matrix 

Since we expect to find two static-light mesons at large R in our Gmm correlator, we look for a positive eigenvalue 
"SS" spectral component corresponding to two S'-wave mesons and a negative eigenvalue SP component corresponding 
to an 5-wave and P-wave meson. For most of our analysis we omit the PP component, a choice that we justify as 
follows: Squaring the static light propagator suggests that the PP channel would contribute a second nonoscillating 
spectral component at a higher energy (about 0.5 in lattice units) than SS and with a smaller amplitude (about 0.2). 
We present results, however, that include a excited state component inspired by the PP contribution. The SS and 
SP components would be expected to have a smooth dependence on R for large R. To these two spectral components 
we add a third, corresponding to a conventional Wilson-loop contribution at short distance R. With the staggered 
fermion phases included, the net Wilson-loop phase factor (— l)(-f^+i)(^+i) produces a transfer matrix eigenvalue with 
a phase {—)^^^ that oscillates with R. 

Our proposed fitting ansatz is thus Eq (|l|) with = 3 and with the explicit SS, SP, and flux-tube eigenvalues 
(respectively) 

Ai(i?) = e-^i(«) 

\2{R) = -e-^^(^) (2) 
A3(i?) - (-)^+ie-^3(fl) _ 

Other components can be readily included. With our choice of sources and sinks the correlation matrix is found to be 
real, so we may take real Z factors. An ambiguity permits changing the sign simultaneously in Zpi and Zmi, which 
we resolve arbitrarily by requiring Zmi to be positive. At large R we expect Vi{R) to approach 2aEs and V2{R) to 
approach aEs + o.Ep and at small R we expect V3 (i?) to correspond roughly to a Coulomb plus linear heavy quark 
potential. At intermediate R we expect mixing among these states. Avoided level crossing may occur at even R 
between intermediate states 2 and 3 and at odd R between intermediate states 1 and 3. 
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IV. RESULTS 



A. Potential 



With the factorization inherent in our ansatz and the choice TV = 3, we are fitting three correlators with nine 
parameters. Our fitting range in T varies over the data set as shown in Table | with typically f or more degrees of 
freedom. The goodness of fit supports our ansatz. 

Our selection of fit ranges compromised between the need to get acceptable fits and our intention to vary the end 
points trnin and tmax smoothly as a function of R. At low R the Gff (flux-tube-type) correlator has quite small 
errors, whereas at larger _R, errors increase, especially at higher T. Thus for small R we can set a higher t,„in in the 
Gff correlator to improve chi square without loss of information. At larger R there seems to be little improvement 
in goodness of fit in setting a high ^min, and the large statistical errors at high T give no advantage to setting a higher 

^max- 

To give an impression of the quality of the fits, we plot the absolute value of the correlation matrix elements vs T 
for two values of R in Figs. ^ and ||. Also plotted are the absolute values of the fitting functions. The errors on the 
observed Gfm and Gmm correlators also give an indication of the signal obtainable with the random source method. 
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FIG. 4. Absolute value of the correlation matrix elements vs t at 7? = 3. The lines connect the best fit values. The fit ranges 
are [4,9] for Gff, [1,9] for Gfm and [2,9] for Gmm- 
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FIG. 5. Same as Fig. |, but with R = &. The fit ranges are [2, 9] for Gff, [1, 7] for Gfm, and [2, 7] for Ga 
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Results at i? = 3 are a good representative of the small R correlators and the degree to which the fit results are 
affected by the choice of end points. Fitting the three correlators Gff, Gfm, and Gmm over the ranges T € [6,9], 
[2,9], and [1,9], respectively, gave Vi = 1.60(5), V2 = 1.89(11), and V3 = 1.026(6) with x'^/df = 14.6/12. Changing 
the fit ranges to t e [2,9], [2,9], and [2,9] increased x^/df to 23.6/15 and gave Vi = 1.62(5), V2 = 1.88(11), and 
V3 = 1.025(8). 

To see the effect upon the mixing analysis of including other states, we experimented with adding an excited state 
modeled after the two- meson PP spectral component, which contributes in the same way as the SS component. The 
fourth spectral component is denoted 14(7?). Doing so increases the parameter count to 12. To assure stability of 
the fits, we fixed the two-meson energies Vi{R) — 2Es and V4{R) = 2Ep, leaving 10 free parameters. We found 
acceptable fits. The unconstrained energies V2 and V3 agreed within errors with results from the three-spectral- 
component ansatz. For example at i? = 3 over fit ranges t E [4,9], [1,9], [1,9] we find x^/df = 22.9/14 with 
V2 ^ 1.95(2) and V3 = 1.024(4). 

Our fit results for the three-spectral-component ansatz are listed in Tables || and ^ and plotted in Fig. |^. For 
small distances the flux tube energy V3{R) is smallest and the flux tube state dominates the large T behavior of the 
correlation matrix, while at large distances the two- meson energies Vi (i?) and V2 (R) are smaller and the two- meson 
states dominate the correlation matrix at large T. With our choice of light quark mass the first level crossing occurs 
at i? = 6a, or 0.98 fm. The string is broken. It is interesting that the energies Vi{R) and V2{R) are very nearly equal 
to their asymptotic values throughout. Thus we see no spectral evidence of a meson- meson interaction at the level of 
our statistics. 

It is clear from these results that mixing between the flux-tube and two-meson channels is weak. There is no evident 
rounding of the potentials normally associated with avoided level crossing. At higher order in mixing we would expect 
to require two- meson spectral components in the flux tube Gpp correlator. They should appear as a result of the 
breaking and rejoining of the string. However, the amplitudes for both terms in this correlator are small enough that, 
if it weren't for the enforcement of a common spectrum and factorization in our fit ansatz, they might have been 
missed. The converse presence of the string term in the diagonal two-meson correlator Gmm can be accounted for by 
the "box" diagram in the quark correlator that resembles a Wilson loop. 

As a check of mixing between the flux-tube level and the two-meson levels, we examined the transition amplitude 
Gfm to see if, by itself, it contained both types of spectral components [|2^. To do so we carried out a separate 
three-exponential fit to the transition amplitude Gfm alone, fixing the energies of the SS, SP, and fiux-tube spectral 
components to the values found in the multichannel channel analysis, but adjusting their amplitudes for a best fit. 
For r < 5, where the spectral components are clearly nondegenerate, we found that the amplitudes for the SS and 
fiux-tube components were both nonzero at the three to five-sigma level. Thus our results, combined with unitarity, 
confirm mixing and imply string breaking. 

We have required that the valence and sea quarks match, thus assuring that the correlation matrix is a power of the 
transfer matrix. However, if mixing effects are very weak, had we chosen, instead, to omit the dynamical sea quarks 
altogether, it would very likely be difficult to detect the consequent inconsistencies. We analyze mixing further in the 
following subsection. 
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FIG. 6. Heavy quark potential and first two excited states vs separation R. The dashed and solid lines give the asymptotic 
values 2aEs and Ep + Es- Jackknife errors are shown. 



B. Modeling the mixing 



Drummond, Pennanen and Michael analyze their results in terms of a transfer matrix model that mixes the two- 
meson and flux-tube states |0,|l6|. Our approach differs slightly, because our multi-exponential fit carries more 
information and because our multi-channel contributions complicate the analysis. A perturbative model of string 
breaking starts with a zeroth order pure flux-tube state and a pure two-meson state, consisting of either a pair of 
unperturbed S'-wave mesons (S'^S^) or unperturbed S- and P-wave mesons (S^P^). (By extension, we could include 
the P^P^ channel.) At a given R the lattice mixing between the states can be described by a transfer matrix on the 
unperturbed basis 

( \\{R) X \ 
T{R) = \1{R) y (3) 
V ^ y >S,{R)j 

that evolves a state across a Euclidean time slice. The rows and columns are arranged in the order S'"S'", S'"P°, and 
F'^ . We assume a small value for the mixing x between P° and S^S^ and y between F° and 5""^". Although the 
S^S'^ and S^P^ states may mix, for simplicity, we have ignored this effect. The diagonal elements correspond to our 
conventions for our fit ansatz (|^). The 2x2 correlation matrix connecting our flux-tube and two-meson source and 
sink states at a given R is 

G(i?, T) = Z"{R)T{Rf+^Z'^{R) , (4) 

where is the unperturbed 2x3 matrix used in our flt ansatz. The diagonal elements of G(T) and the potentials 
are unperturbed at flrst order in x and y. The off-diagonal correlator to flrst order is 

GFM{R,T)=A,[X',iR)f+'+A2[X'^{R)f+'+A,[Xl{R)f+' , (5) 
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where 



Ai — Zfi{R)Zmi{R) — Zp^{R)Zlj^{R) 



xuR)-xm 



Ai - Zfi{R)Zm2{R) - Z%j^{R)Z\i2{R) yQ{ji\^ _ \^(B) 

A3 = Zf3{R)Zm3{R) = -{Ai+ A2) 

To first order in the mixing parameters, the Z factors are unperturbed, and we may equate Z — Z'^ to obtain three 
constraints for the two mixing parameters x and y. The third constraint is a sum rule. To apply the model we chose 
to impose the sum rule as an a posteriori constraint on the parameters of the fit for each R: 

f{Z) = ZpiZmi + Zp2Zm2 + Zp^Zm3 = , (7) 

and use the first two conditions to determine the mixing parameters: 

x=[X\{R)-\l{R)]ZFilZF:i 

y=[\l{R)-\l(R)]ZF2lZF3 . (8) 

This was done by linearizing the sum rule in the vicinity of the minimum of the unconstrained x^- 

/(Zo) = = /(Z,) + V/(Z,) • (Zo - Z,) . (9) 

It is straightforward to determine the attendant increase in , the shift in parameters, and the decrease in errors. 
Since this procedure assumes the sum rule can be linearized, it is valid only to the extent that the increase in is 
small. 



Table III lists results. Shown are the values of the coefficient ratio 

d = —{ZfiZmi + Zf2Zm2)/ZfzZmz , (10) 

before and after imposing the linearized sum rule constraint as well as the shifts in the values of x and y 

obtained after imposing the constraint. The ratio d should be one if the sum rule is satisfied. We see that the 
increase in is smallest for R > 5.66, but it is otherwise unacceptably large. Thus the mixing model suits our 
three-exponential ansatz only for larger R. The agreement improves considerably when we include the two-meson 
PP spectral component, as discussed above, and fix the two-meson energies Vi{R) = 2Es and V4{R) = 2Ep, leaving 
10 free parameters. Although we include the fourth component in the fit, we still consider only a three-state mixing 
model. (In effect, we have set mixing to the fourth level to zero.) Results are shown in Table Now the mixing 
model seems plausible for i? > 2. 

The mixing model makes separate predictions for the connected and disconnected meson-to-meson correlators, 
which provides an additional constraint on the mixing parameters. As a test of the systematic error arising from 
model assumptions, we have tried refitting all of our data to a purely four-component mixing model, with separate 



disconnected and connected correlators. While the resulting mixing coefficients repeat the trends of Tables |I| and 



[V, the values differ as much from those of the tables as the two tables do from each other. 
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FIG. 7. Absolute value of the mixing parameters x and y vs separation R. Odd and even series are distinguished. 

If we now take the mixing model at face value, it is interesting to consider how the strength of mixing varies with 
R. In Fig. we see a significant decrease in both x and y with increasing R, which is to be expected partly from 
as the eigenvalues cross. There is also a pronounced difference in mixing strengths at even and odd R, suggesting 
a suppression of mixing between unperturbed oscillating and nonoscillating levels. Obviously, given the coarseness 
of our lattice, the comparison with a model Hamiltonian must be done judiciously. However, even these first, crude 
QCD results should help constrain the phenomenological analysis of quarkonium decay p2[. 



V. CONCLUSIONS 



We have studied string breaking in the heavy quark interaction at zero temperature in QCD. Our calculation used 
two flavors of light staggered quarks and gauge configurations generated in the presence of the same quarks. We 
extended the analysis of the spectrum of the transfer matrix in the staggered fermion formalism to treat our nonlocal 
sources. By adding explicit two-meson states to the conventional flux-tube state, we obtain the expected result that 
the two-meson state is energetically favored at large distance. Our two-channel correlators fit a model with three 
factorizing spectral components and a partially constrained, extended model with four spectral components. Our 
results are also consistent with a simple three- and four-state model with weak mixing. The mixing coefficients in the 
mixing model matrix appear to decrease at the level crossing points. Our principal finding is that mixing between the 
flux-tube and two-meson channels is indeed weak. Thus we see why string breaking has been missed in the flux-tube 
channel by itself. While our matching of dynamical and valence quarks is designed to satisfy unitarity, with our 
statistics we have not found compelling evidence for quark loop effects. Doubtless, our results could have just as well 
been reproduced in a quenched simulation, where inconsistencies with unitarity would appear only at higher order 
in the mixing matrix elements. However, we also find that the transition correlator by itself connecting fiux-tube 
and two-meson channels is nonzero for all r and shows both string-like and two-meson-like spectral components, at 
least for r < 6. Thus, our results, combined with unitarity, require string breaking. It is recommended that future 
staggered fermion studies closer to the continuum take care to restrict the two-meson wave function to the light-flavor 
singlet channel. 
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APPENDIX A: RANDOM SOURCE ESTIMATORS APPLIED TO STRING BREAKING 



We use random source methods to calculate the all-to-all quark propagators in this study. Here we outline the 
method as applied to string breaking. We begin with the conventional light quark propagator 

Gf{x,y) = {q''{x)^{y)) , (Al) 

which satisfies the equation 

(^ + ml)°f,G^^(2/,^) = ^-^,,, . (A2) 

Summation over repeated indices is implied unless noted otherwise, p denotes the usual staggered lattice Dirac 
operator. 

The static (infinitely heavy) quark propagator can be found from an expansion to leading order in 1/m. The result 
depends on whether propagation is forward or backward in time. For t> t' , omitting the heavy quark mass factors, 

Gt[{x, t), {y, t')] = {h%x, t)h\y, t')) = 6s,y ^ f[ ul{x, r) j , (A3) 

where 

Ulix, t) = a{x, r; fJ,)Ul{x, r) (A4) 
and the a{x,T;ij) are staggered fermion phase factors. For t <t' 

(t' 1 \ ab 

l[UMr)\ , (A5) 

where 

Uij.{x, r) = -a{x, r; /u)i7^(f, r) . (A6) 
We will be interested in static-light mesons, which we construct from the Grassman fields as 

B{x, t) = Y,Y. ^)'^"'(^' y^i'^y^ - 

V o,,b 

B{x, i) = E E *^pt\y^ ^)^'(^' • 

The wave function pt is hermitian, pj = pt, and depends via gauge fields on t. We use 

3 ( /ro-l /ro \ «^ ] 

pf{x, y) = 25'''5s,y + + Ijl, t) j 5s,y-r,[. + ^l^^ ' ^Z^' *) j -^-.ff+^oA | • (A8) 

The space-like gauge-link matrices can be taken as (APE) smeared gauge fields. Thus the tilde here represents both 
smearing and the inclusion of the fermion hopping phases. On the other hand the time-like gauge-link matrices are 
not smeared. 
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For the static-light meson correlation function with source at the origin we obtain 

GeiT) = (i?(0,r)S(0,0)) = Gf [(0,0),(0,T)]p^^(0,f)Gf [(f,r),(y,0)]pg''(y,0) 

= W^i6,0)V^'[i6,0),i6,T)]p'^^i6,x)Gfm^^^ . (A9) 

where 

/t-i \«'' 

V^'lix, 0), (f , T)] = ill '^^l^- r)]\ (AlO) 

is the product of link matrices and hopping phase factors from {x, 0) to (x, T). In the last line we introduced Gaussian 
random numbers, {'ri^"'{x)rf'{y))n — 5°'^5x,y, to compute the trace. Multiplied by pt they form the source for the light 
quark inversion and are then used again in the computation of the static-light meson correlation function as seen in 
eq. ( |A9| ). Introducing the smeared, random source propagator 

= pl\x,y)G\%y,t), {z,t'Mf{z,u)r^\u,t') , (All) 

we can write the static-light meson correlation function as 

GsiT) = (7,ta(o,o)F'^''[(0,0),(0,r)]*f,(0,T)), . (A12) 

In numerical practice this result is computed for a source at any lattice location and averaged over the space-time 
volume. 

For the computation of the heavy quark potential and the investigation of string breaking we will need a heavy 
quark-antiquark "string state" 

OF{S,y,t)=h''{x,t)V''''[{x,t),{y,t)]h''{y,t)/^,, (A13) 

where V"'''[{x,t), {i/,t)] is a superposition of products of (APE smeared) gauge fields in time slice t starting at x and 
ending at y. 

We are then interested in the diagonal correlator between a string state at time and one at time T —we consider 
the connected part only 

{Of{R, 0, T)Of{0, R, 0)) = W{R, T) , (A14) 

with smeared space-like gauge field product segments in the time-like Wilson loop W{R,T). Except for the fermion 
hopping phases, which give a net factor (— )(^+i)('^+i) independent of staggered fermion phase convention, this is 
the correlation function usually considered for the computation of the heavy quark potential. As usual, in practice 
this quantity is averaged over all choices of lattice origins and on-axis displacements R. We also computed it for two 
off-axis displacements, namely R — \/8 and 2^/8. For off-axis displacements V°'''[{x,t), {y,t)] is constructed from a 
symmetric set of space-like paths joining the endpoints. 

We are also interested in the off-diagonal correlator between the string state and the two-meson state 

Oa/(0, R, t) = B(i?, t)B{Q, t) . (A15) 



{Of{R,0,T)Om{0,R,0)) 
= V-'[{x, 0), {x, TW'HR, T), (0, r)]F^'^[(^, T), {R, 0)]pg^(0, x)Glf[{x, 0), {y, 0)]p^"(y, R) 
= (v^^iR, 0)V'''[iR, 0), (i?, T), (0, T), (0, 0)1*^,(0, 0)), . (A16) 

The V with four coordinate arguments stands for the product of time-like and (APE smeared) space-like gauge fields 
that make up the part of a time-like Wilson loop without the spatial segment at time t = 0. Again, this correlator 
is computed with the Gaussian random source method, in a manner similar to the static-light meson correlation 



function (A£). The other off-diagonal correlator, {B{R,T)B{0,T)S{0, R,0)) , is the hermitian conjugate of ( [AT 
This quantity is averaged over the same choices of origin and displacement as the Wilson loop. 
Finally, we need the two-meson correlator 
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{Oiji6,R,T)OMiO,R,0)) 
= V<^\R, 0), (i?, T)]p\'{R, S)Gf[{x, T), (y, T)]p^^(y, 0) 
X T/^/[(0, T), (0, OM^O, z)Gf[{z, 0), (.1, 0)]pg"(", R) 
- ^'^''[(i?, 0), (i?, T)]pif{R, x)Gf{ix, T), (y, 0))pf(''(y, R) 
X F^/[(0,r),(0,0)]p^^(0,z)Gf ((z,0),(u,T))p^^(u,0) . 



(A17) 



We compute the two terms in eq. ( A17 ) again with the Gaussian random source method using two independent 
Gaussian sources, rj and C, one for each of the two hght quark propagators: 

{Oli6,R,T)OMi.0,R,0)) 
= (*^(0, oW'^iR, OW'^'iiR, 0), {R, nU'fliR, r)Ct^(0, T)V'f[{0, T), (0, 0)])c 
- (77^'^(i?, 0)V^\{R, 0), (i?, T)]^\{R, T)),(Ct^(0, r)l/'=/[(0, T), (0, 0)]vI/{(0, 0))c • 

Again, this quantity is averaged over the same choices of origin and displacement as the Wilson loop. It is important 
to note from our expressions for the correlators that all can be computed using 0{N) methods for N random sources. 

APPENDIX B: TRANSFER MATRIX APPLIED TO STRING BREAKING 

We review the Sharatchandra-Thun-Weisz (STW) Fock-space formulation of the staggered fermion partition func- 
tion with specific application to the operators used in our string breaking study . 



1. Introduction 

The fermion action is given by 

S = ^5*3(9*, gt,[/t) +^ar,t;4 qr,t+lUr,tqr,t- qr,tUl~qr,t+l 
t r,t 

where the spatial part of the action is 

Saiqt^qt, t^t) = X! 1 '^'^qr.tqr.t + ^ ar,t;i qr+i,tUr,t;iqr,t " qr,tUlt;iqr+i.t 
r \, i 

With the STW phase convention 



(Bl) 



(B2) 



(B3) 



With this choice the phases have no t dependence. The imaginary time variable ranges over 0, . . . , — 1, and the 
antiperiodic boundary condition requires qr.N — ~<7r.o- 
The partition function is given by 



^~ J dqdqexp{S) 



(B4) 



where q and q denote the full set qr,t and gr,t- We would like to convert the Grassmann integral into a Fock-space 
operator trace. To this end STW first eliminate the KS phase in the time direction altogether, by changing variables 
to = ar^t;4q- The spatial action becomes 



Ssiqhqt, Ut) = ^ \^mqltqr,tar,t-i + ^ a'r^ui ql+i,tUr,t-rqr,t ^ ql^tUl^t-tqr+i,t I 



(B5) 



where 
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(B6) 

Then STW introduce a dummy set of Grassmann variables pt = ql and p] — qt- The action is then 

S^Y. S3{qlpl)/2 + S3{pt,qt)/2 - [plpt+i + qUt+l] ■ (B7) 
t 

where we have suppressed the sum over spatial coordinate r. The partition function now includes integration over 
the dummy variables with a delta- function constraint: 

Z = J dqUqdpUpS{p-q^)S{q~p1)exp{S) . (B8) 

Here the delta function implies a product over delta functions on each lattice site. The Grassmann delta function is 
simply 

S{q-p^)^iq-p^) . (B9) 

We introduce a Fock space by associating creation and annihilation operators with each Grassmann variable. The 
corresponding Fock space operators are denoted by a hat. It is convenient to introduce Grassmann coherent states 
11: 

q\q) = gk) p\p) =p\p) 

(gt|^t ^ ^^t|^t (pt|pt ^ ^pt|pt , (BIO) 
The coherent states satisfy completeness and trace relations: 

1 = y dqUq \q) {q\ expi-q'^q) 
TtA^ J dqUq{q'<\A\q)exp{q'<q) . (Bll) 

With these identities we can reorganize the factors in exp(5) in the form 
exp{S)6{p~q^)6{q-p^) 

= [| {e-f^f*+i-«^«'+i ^qIpI e^3(«^'^^*)/2 : (p_qt)(^_pt) . e^'^^'^'P^^^ \qtPt)} ■ 
t 

where :: denotes operator normal ordering. Using the completeness and trace identities and the antiperiodic boundary 
condition, we can then write the partition function in terms of a transfer matrix operator: 

Z = TrT^ , (B12) 

where 

r = e^^(«''^^''/2 : (gt -p){q- pt) ; e^3(?.p)/2 (2^3) 
This is a manifestly Hermitian, but not positive definite, transfer matrix. 

2. Quantum mechanics 

To see how this transfer matrix works, as a warmup exercise, suppose there is only one spatial site (quantum 
mechanics). The transfer matrix is then 

r = (1 + mgV) . (^t - p)(g - p) : (l _ mqp) . (B14) 

It has eigenstates (m = sinh7) 

|0) , |o) , |±c) = (1 T e^Tgtpt) |0) /^l + 6^27 (B15) 
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with eigenvalues 1,-1, zte^^. Thus the ground state is |— c). The "degenerate" p and q states arc interpreted as 
single particle and single antiparticle states, and the state |+c), as a state with a particle and antiparticle pair. The 
partition function is {N even) 

Z = 2 + 2 cosh(A^7) . (B16) 
Next we consider the propagator for a single massive fermion, created at time slice and propagating to time slice 

T. 

G{T) = j dqUqeMS)qTql/Z . (B17) 

Converting to the Fock space basis as we did with the partition function leads to 

G{T) = Tr[T^-^-^p^T^+^p]/Z . (B18) 

This example shows that a T + 1 power is natural. If we assume that N is large, so only the ground state contributes 
to the traces, we have (T <C N) 

G{T) = e-T'^/[2 cosh(7)] (B19) 

with 7 interpreted as the energy of the propagating state. For large mass the correlator is approximately l/(2m)(-^+^^. 
The antiparticle propagator (with q^^qo instead) is also easily computed with the result 

G{T) = -(-)^e-T^/[2cosh(7)] . (B20) 

Next, we introduce a second heavy flavor of mass M {M = sinhF), denoted by h. We consider the propagation of 
a heavy-light meson. For simplicity we use an interpolating operator with a local (point) wave function. In this case 
the Grassmann integration involves the combination qj.hTh\qQ. Following the same steps as before leads to 

Gb{T) = -(-)^e-(T+'^)^/[4cosh(7) cosh(r)] . (B21) 



3. String-breaking operators 

Here we construct the operators required for our string-breaking study. Two interpolating operators, 0\ and Oj^^, 
are used, one that creates a static quark- antiquark pair with a connecting flux tube and another that creates a pair 
of static-light mesons. A simple way to construct the flux tube operator is to introduce a new heavy flavor H of 
mass Mh- Our conventions for the on- axis (direction i) string-creation operator can be obtained from the Grassmann 
product {2MH)^ho,tHo,tHR^thR,t after integrating out the new flavor in the static limit Mh oo. The result is a 
heavy quark- antiquark creation operator connected by the static quark propagator. In the Grassmann basis it is 



Oy{R,t)=F^R,t) = ho,, 



R-l 



r=0 



hm,t 



(B22) 



Our corresponding annihilation operator is similarly generated from the product {2MH)^hii^tHR^tHo,tho,t, yielding 



OF{R,t)=F{R,t) = hm,t 



r=0 



(B23) 



(With all of our static quark propagators, it is convenient to introduce a separate, heavy flavor for each straight-line 
segment.) The meson-antimeson creation and annihilation operators are simply 

Oi^(ii,t) = B(0,0)B(i?,0) = hQ,tqo,tqRt,thRi,t 
OM{R,t) = B{R,t)B{0,t) = hRi,tqRi,tqo,tho,t ■ 



We are interested in the two-channel correlators in the static limit Mh 



(B24) 

oo with the heavy quarks fixed at and R: 
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GAsiR, T) - {Oa{R, T)OliR, 0)) {2Mhf'^ , (B25) 



where A,Bg {M,F}. The conversion to Fock space proceeds as before. With 6{g — h^)6{h — g^) the Fock space 
operators are 



R-l 



n {-ar,t■,^Ult■,^) 



r=0 
R-l 



l'Ri,t- 



t-A 



W {ar,t;iUr,t;-, 



r=0 



(B26) 



OliiR) = go.taa.tAqo,tPRi,taR.i,t;4:f^Ri,t 
OMiR) = hm,t(^Ri,t-APmAo,t(^o,t-49Q.,t ■ 

We recah that Q!r,t;4 is independent of and ct^^-^ = (^)'''; so the dagger indeed denotes the Fock space 

hermitian conjugate. The desired correlators are then 

GAB(i?,T) = (vac|0A(i?)T^+'0l3(i?)|vac) . (B27) 
An eigenstate \n) of the transfer matrix with eigenvalue A„ contributes 

GABniR,T) = Z*A^iR)ZBniR)[XniR)f+^ , (B28) 

where 

ZAniR) = {n, R\ OAiR) \vac) . (B29) 

This result is the basis for Eq. (|^). 

4. Staggered spin and flavor considerations 

The external states in our analysis are built from two operators: one that creates a static quark-antiquark pair at 
separation R and one that creates a pair of static quark-light quark mesons at separation R. Both the static quark 
and light quark carry four continuum flavors. One may ask, in the continuum limit, what spin, parity, and flavor 
combinations occur? 

a. Flavors and spins of the static-light meson 

Golterman and Smit [p5| give an analysis of the flavor content of staggered quark-antiquark mesons. Our static-light 
meson wavefunction has a local component (zero displacement) and a component with displacement 2 along any axis. 
An even displacement does not change the light-quark flavor and spin wavefunction, so for simplicity we analyze only 
the local component, 

h{2y + rf,t)q{2y + rf,t) , (B30) 

where the cm. offset ff ranges over the eight sites in a unit cube. 

The more familiar flavor-spin content is displayed in the notation of Kluberg-Stern et al. : 

g(2y + 77,i) = ir^V'^(y) (B31) 
where a sum over four flavor a and four spin a indices is assumed, and 

^T = ^TiTiTiT- (B32) 

A similar expression holds for the heavy quark operator h{2y + rj^t). We have ignored the gauge connection in these 
expressions. The coefficients F^" give the flavor-spin wavefunction of a quark created at site rj. 
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The zero-three-momentum projection of the local static-light operator is given by 

Bfj{t) = J2H'^y + V,t)qi2y + ri,t) . (B33) 
y 

All eight operators, corresponding to the eight values of ff, belong to Golterman's class 0. Because the static quark 
must propagate in place, the correlation matrix for these eight states is diagonal in and independent of 77. Apart from 
a volume factor, it is exactly the static-light correlator we have computed. Linear combinations of these operators with 
coefficients differing only by signs yield operators belonging to the rest-frame symmetry group of the discrete transfer 
matrix. For example, the operator belonging to the one-dimensional representation 1"' in the Golterman-Smit 
notation is given by 

In the Kluberg-Stern et al. basis this operator is written as 

^ ^(2^) X 75 ± 70 X jQ]q{y) (B35) 
y 

where the first gamma factor in the tensor product operates on the spin basis and the second, the flavor basis, and 
the sign is plus for even t and minus for odd. 

This operator generates both a pion-like and sigma-like state in the continuum, in which, respectively, the light 
quark is found in an S or P orbital around the static quark. Similar linear combinations yield eight operators, 
altogether belonging to the irreps 1^^, 1^ , 3""++, and 3""++, each interpolating one of eight S'-wave channels: 
two pion-like (_B-mesons) and two rho-like (i?*-mesons). Each of these is paired with one of eight P-wave channels: 
two sigma-like and two pseudovector-like, respectively. All S'-wave states are degenerate, as are all P-wave states. A 
fixed offset ff corresponds to a distinct linear combination of these degenerate states. Alternatively, we may say that 
a fixed offset ff corresponds to the pairing of a single flavor-spin species of light quark and heavy quark. 

Thus our two-meson operator interpolates all combinations of B and B* and their P-wave counterparts. Our 
procedure for computing correlators sums over all cubic rotations of R at fixed |P|, so is consistent with zero total 
angular momentum for the two-meson and flux-tube states. Transitions of the type SS — > SP are accompanied by a 
change in orbital angular momentum. 

Static-light mesons in the Golterman classes 1,2, and 3 are similarly degenerate within each class with multiplicities 
24, 24, and 8, respectively. One would expect in analogy with the light qq mesons that on a coarse lattice the energies of 
the class nonzero static-light mesons are higher than that of class 0. In the continuum limit all 64 become degenerate 
because of flavor and heavy quark symmetry. A multiplicity factor of 16 comes from separate heavy-flavor and 
light-flavor symmetry and a multiplicty factor of 4 = 1 + 3 for degenerate singlet and triplet spin combinations. 



b. Flavors and spins of the static-static meson 



The flux-tube operator is designed so that it belongs to the same representation of the symmetry group of the transfer 
matrix as the meson-antimeson operator, thus permitting mixing of the states they create. The zero momentum 
projection of the mcson-antimeson operator at any separation R is in Golterman class 0. At P = the flux-tube 
operator creating a static quark-antiquark pair is trivially one of eight Golterman class operators. For it to remain 
in this class as the quarks separate, one must include the hopping phases ±a^{x) with the gauge- link operators that 
excite the electric flux accompanying the creation of the static pair ( A4 ) . A simple way to see this is to consider that 
the operator class is a symmetry of the transfer matrix, so acting upon the zero separation state with the heavy-quark 
hopping matrix, which incorporates the phases, preserves the operator class. The eight operators consist of two pion- 
like and six rho-like operators. More specifically, the static-quark-antiquark operator at all P is a linear combination 
of the two local (heavy quark) pseudoscalars tts and ttqs and the six local vector mesons pi and pio- In the 6-quark 
system, these states would be classified as the T and the rjb- 
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In free (static) propagation the even parity P-wave partners of both are trivially absent in the static limit. On 
the other hand, in the presence of dynamical staggered quarks and in the strong-coupling limit, we see no reason 
that dynamical quark loops could not generate an anomalous, opposite-parity partner in the Wilson loop correlator at 
second order in the mixing coefficients. Thus, the appearance of both oscillating and nonoscillating components in the 
Wilson loop would be a signal of string-breaking at strong coupling. In the continuum limit the anomalous components 
should vanish, in keeping with the expectation that lattice staggered and nonstaggered fermion actions have identical 
limits. On our moderately coarse lattices, in the transition correlator connecting the flux-tube and meson- antimeson 
channels, the signal appears in leading order, but from our mixing analysis, we find that the anomalous parity 
component is smaller. 



c. Light flavor artifacts in the meson-meson correlator 



The peculiarities of staggered flavor symmetry require paying special attention to light quark flavor counting. Two 
issues confront us. First, in our meson-antimeson channel, in the continuum limit the light quarks can combine as a 
flavor singlet or flavor non-singlet with, in principle, distinct energies. Since the flux-tube state is necessarily a light- 
flavor singlet, only the singlet components mix, leaving the nonsinglet as a potential additional spectral component. 
Second, the gauge configurations were generated in the presence of Nf ~ 2 flavors, through the usual device of taking 
the square root of the intrinsically four-flavor fermion determinant. However, the light quarks in our source and 
sink carry four flavors. A mismatch between the number of valence and sea quarks gives rise to additional spectral 
components in the disconnected meson-to- meson correlator. 



To count flavors, we construct a toy model similar to the mixing model of subsection [VE. However, for present 
purposes we ignore negative transfer matrix eigenvalues and deal directly with the Laplace transform of the correlators. 
Further, it is sufficient to concentrate on the meson-meson correlator Gmm- 

We begin with the disconnected meson-meson correlator Gmmd at a fixed separation R. As we have noted above, 
the static- light mesons are created in Golterman class 0. In the absence of glueball exchange between the meson 
and antimeson, they remain in class 0. Glueball exchange can excite class n > levels in pairs. Such an effect is 
presumably short range, because the glueball mass is of order one GeV. Thus, on a coarse lattice we expect four 
distinct energy levels, and the Laplace transform of the Gmmd correlator has the form 

3 ^2 

Gmmd = V " ■ (B36) 

n— 

At least for R larger than a few tenths of a fermi we expect the lowest level to have the largest overlap with the source 
and sink so that >> Zn for n > 0. So we simplify our model, writing 

I"! 

Gmmd — + . . . . (B37) 

tj — tjQ 

where the remaining terms are higher energy and weakly coupled, except, possibly, at small R. 

In the continuum limit flavor symmetry requires only two energy levels, depending on whether the light quarks form 
a flavor singlet or nonsinglet. Thus we use to denote the singlet and the nonsinglet. Since our source/sink 
operator generates a single flavor-spin species of light quark and light antiquark, for four light flavors the relative 
weight of singlet to nonsinglet is 1 : 3, and we get 

Gmmd Icont = -p, + "^1? 5^- (B38) 

Next, we consider the connected correlator Gmmc- We treat light quark pair annihilation and creation as a weak 
process with amplitude x and analyze this correlator as a perturbation series in x. Thus the leading term in the 
connected diagram has the form 

^MMC - _ £;^)2(^ _ E^.) ^^■^^> 

where Ef is the energy of the flux-tube state. The amplitude x gives the weight for annihilation or creation of a 
single light quark/antiquark flavor. Here is the sequence of events: The source state couples to the lowest level, which 
propagates with energy Eq. It then annihilates with amplitude x to produce the flux-tube state, propagating with 
energy Ej. Pair creation leads back to the lowest level, which propagates with energy Eq before coupling to the sink. 
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The next higher order contribution to Gmmc comes from a quark-loop insertion. The sequence of events is the 
same as in the leading-order contribution, except that an extra pair creation and annihilation process occurs in the 
flux-tube state. Pair creation on a coarse lattice leads to any of the four levels Ei at the same order in x. In the 
continuum limit and large R these levels are degenerate, and would each count with weight x'^Nf/A^ the factor -/V//4 
arising from the fcrmion determinant. For simplicity, to model the full effect of the internal quark loop, we assign it 
a weight x'^ with x'"^ = x'^Nf in the continuum and represent only the lowest level Eq. Thus we have 



^MMC - ^MMC {^E~Eo){E- Ef] ■ ^ ' 



Continuing with these simplifications to all orders, we may sum the perturbation series to obtain Gmm = Gmmd + 
Gmmc-, where 

^ ^■ A + {x^-x'^)/[{E-E,){E-Ef)] 

^^^ = ^ E-E,-x'y{E^Ef) + ■■■■ ^^^^^ 

The ellipsis represents contributions to the disconnected diagram from light-flavor nonsinglet terms. 

This toy flavor-counting model shows the desired shifted pole at Eq + x''^ /{E — Ef), but additional, comphcating 
spectral components at Eq, Ef, and the light-flavor nonsinglet level(s). These additional components are an artifact 
of our choice of source and sink. While, in principle, one may carry through the spectral analysis of the correlators, 
remembering to include the artifacts, the analysis would be unnecessarily complicated. For the present numerical 
study, weak mixing and strong flavor symmetry breaking on our coarse lattice render the extra spectral components 
harmless. First, strong coupling breaks flavor symmetry and lifts the levels Ei, . . . above Eq, making it less likely 
to confuse them with Eq. Second, because of weak mixing, we are unable to detect the internal quark loop directly. 
Thus we measure the mixing parameter x, but do not see x' . 

Closer to the continuum, it is desirable to replace the meson-meson interpolating operator with a better one, building 
in an explicit projection of the light-quark flavor singlet. An example is a state formed from the tensor product of 
a a meson and a static quark-antiquark flux-tube. Clearly, such a projection eliminates the nonsinglet contribution. 
But to eliminate the other artifacts also requires x^ = x'^ in the toy model, which is not automatic when there are 
two sea quark flavors but four source and sink flavors. The remaining artifacts are eliminated by adjusting the weight 
of the disconnected diagram, relative to the connected diagram. To be precise, once a projection to an SU(4) flavor 
singlet has been done in the state M, the proper weighting is Gmmd + Nf /AGmmc- 
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TABLE L Fit ranges in t, chi square, potentials, and k vs R. Jackknife errors are given. 
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[1,12] 


12 121 


20 0/23 


1.57(3) 




1.71(10) 


5229 1'4^ 


2.0 


[4,11] 


[1,10] 


19 101 


24.0/18 


1.62(5) 




1.91(13) 


81 81 i^ 5"! 


2.83 


[3,10] 


[1,9] 


\2 91 


11.1/16 


1.61(3) 




1.81(8) 


1 nnfi6n 3~i 


3.0 


[4,9] 


[1,9] 


[2 91 


23.4/14 


1.63(5) 




1.91(13) 


1.024(4) 


4.0 


[4,9] 


[1,9] 


[2 91 


16.2/14 


1.61(6) 




1.85(13) 


1 198(10") 


5.0 


[2,9] 


[1,8] 


[2 81 


18.6/14 


1.63(5) 




1.88(11) 


1 353f3l 


5.66 


[2,9] 


[1,7] 


[2 71 


14.6/12 


1.59(4) 




1.74(9) 


1.458(2) 


6.0 


[2,9] 


[1,7] 


[2 71 


14.8/12 


1.54(4) 




1.69(10) 


1.507(4) 


7.0 


[2,8] 


[2,7] 


[2 71 


9.6/10 


1.69(4) 




1.96(12) 


1 648 f 14) 


8.0 


[2,8] 


[2,6] 


[2 71 


13.3/9 


1.61(4) 




1.81(6) 


1.78(6) 


9.0 


[2,8] 


[2,5] 


[2,6] 


9.6/7 


1.49(3) 




1.54(8) 


1.936(14) 


10.0 


[2,8] 


[2,5] 


[2,6] 




1.58(5) 




1.79(14) 


2.09(2) 


TABLE II. Couplings vs R. 


Odd and even 


R values 


are grouped. (The 


displacements with R 


— \f% and 2\/8 have even 


Cartesian components.) 














R 




Zfi 


Zmi 




Zf2 


Zmi 




Zfz Zms 


1.0 




0.162(3) 


14.0(1.0) 




-0.002(3) 


13(2) 




-1.2733(10) 3.356(3) 


3.0 




0.24(3) 


16(2) 




0.005(4) 


19(5) 




-1.582(14) 4.33(7) 


5.0 




0.23(3) 


16(2) 




0.005(3) 


18(4) 




-1.784(6) 2.9(2) 


7.0 




0.9(8) 


16(5) 




-0.00(2) 


21(5) 




-1.8(4) 8(6) 


9.0 




-0.01(3) 


12.1(7) 




-0.007(10) 


9.2(1.4) 




-2.13(5) 0.2(7) 


2.0 




0.0180(14) 


16(2) 




0.005(7) 


19(5) 




-1.443(5) 1.970(8) 


2.83 




0.0160(7) 


15.4(1.1) 




0.051(4) 


16(2) 




-1.540(3) 1.577(9) 


4.0 




0.0141(12) 


16(2) 




0.056(15) 


17(4) 




-1.70(4) 1.19(3) 


5.66 




0.0087(11) 


14.9(1.3) 




0.06(2) 


14(2) 




-1.780(5) 0.71(7) 


6.0 




0.0067(13) 


13.4(1.2) 




0.12(6) 


12(2) 




-1.890(9) 1.0(3) 


8.0 




-0.003(12) 


15.5(1.3) 




-1(2) 


10(20) 




1(2) 9(15) 


10.0 




-0.014(10) 


15(2) 




0.2(2) 


15(4) 




-2.24(6) 3(3) 
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TABLE III. Results of mixing model analysis based on a three-exponential fit showing the coefficient ratio d! prior to 
imposing the a posteriori sum rule constraint and d after, the increase in , and the mixing parameters x and y vs R. This 
mixing model seems plausible only for R > 5.66. Jackknife errors are given. Values for odd and even R are grouped as in Table 



SX 


J' 

11 




Ci 


X 


y 


i.UU 


u.ozi^o ) 


OZU 




n npi(^(''^A 

U.UOOl^O ) 


n 01 Q{ 


o.uu 


U.ODl^o ) 


900 






n 007/'9^ 


^ nn 
o.uu 




yo 




n nnn/'fi^ 
u.uuui^o J 


n 001 c\(Q\ 


7 nn 




n 






n 0004 rs"! 


9.00 


-1(9) 


0.21 


1.002(10) 


0.00090(15) 


-0.0005(5) 


2.00 


0.14(3) 


1300 


2.4(12) 


0.002(3) 


-0.036(8) 


2.83 


0.43(5) 


230 


0.87(10) 


-0.0044(4) 


-0.0077(11) 


4.00 


0.58(10) 


23 


1.08(10) 


-0.0030(6) 


-0.0075(15) 


5.66 


0.79(9) 


12 


1.00(3) 


-0.0016(2) 


-0.0029(3) 


6.00 


0.85(7) 


5.8 


0.9(2) 


-0.0009(3) 


-0.0029(10) 


8.00 


1.2(4) 


1.6 


1.0(6) 


-0.00(2) 


0.01(15) 


10.00 


0.3(2) 


2.1 


0.9(3) 


0.0003(4) 


-0.0006(5) 



TABLE IV. Same as table III but with a four-exponential fit, constraining Vi(R) 
model now seems plausible for R> 2. Jackknife errors are given. 



2Es and ViiR) = 2Ep. The mixing 



R 



d' 



Ax' 



d 



1.00 


0.52(5) 


140 


1.0(9) 


0.10(14) 


0.007(4) 


3.00 


0.4(2) 
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0.94(2) 


0.051(2) 


0.0089(13) 


5.00 


0.6(4) 


1.8 


0.99(2) 


0.0139(6) 


0.0034(7) 


7.00 


1(2) 


0.00 


1.001(9) 


0.0032(6) 


-0.001(4) 


9.00 


0(20) 


0.032 


0.9(6) 


0.2(15) 


-0.1(7) 


2.00 


0.5(2) 


5.8 


1.06(4) 


-0.037(3) 


-0.016(2) 


2.83 


0.61(15) 


9.6 


1.003(2) 


-0.0200(11) 


-0.0129(3) 


4.00 


0.9(4) 


0.037 


1.0000(2) 


-0.010(2) 


-0.0075(3) 


5.66 


0.3(5) 


4.8 


0.96(4) 


-0.0038(11) 


-0.0025(3) 


6.00 


0.4(6) 


0.83 


1.001(2) 


-0.0011(13) 


-0.0027(3) 


8.00 


1.5(3) 


5.5 


1.1(16) 


0.007(6) 


-0.0008(14) 


10.00 


0.3(2) 


1.6 


1(9) 


0.002(5) 


-0.001(13) 
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